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I. INTRODUCTION AND SUMMARY 


Because of the recent interest in high-precision experiments to test 
gravitational theories, a theoretical framework encompassing a very broad 
and populous class of metric gravitational theories has been developed to 
systematize the comparison between theory and experiment. This framework 
has taken the form of a Parametrized Post-Newtonian (PPN) Formalism 
[Eddington (1922), Robertson (1962), Schiff (1967), Nordtvedt (1968b), and 
Will (1971a)], which treats the post-Newtonian limits of metric theories 
of gravity in terms of a series of metric parameters whose values vary from 
theory to theory. Hand-in-hand with the PPN metric go equations of motion 
describing the response of matter to the metric (the usual "divergence of 
the stress-energy tensor vanishes" for stressed matter, or geodesic motion 
for test bodies). The PPN formalism has been used to study the "classical 
tests" of general relativity, to search for new tests [Nordtvedt (1968a, c; 
1970a; 1971a, b), Will (1971b, d), Thorne, Will, and Ni (1971)], and to 
analyse a wide class of contemporary metric theories of gravity [Nordtvedt 
(1970b), Ni (1972)]. 

In this paper we introduce a revised, and in certain respects, extended 
version of the PPN formalism. We view the original PPN formalism as having 
an empirical emphasis or flavor; each parameter was simply a coefficient of 
a different term in the post-Newtonian metric. In the extended PPN we have 
found it desirable to regroup the metric parameters. The new parameters 
(algebraic combinations of old parameters) are now individually related to 
different physical or conceptual aspects of metric theories. In terms of 
these new parameters, classification of metric theories of gravity becomes 
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simpler and more concise. Every metric theory of gravity can be categorized 

"''Exception: Whitehead's (1922) theory is too complex to fit into the PPN 
formalism. However it has been shown to disagree violently with gravi- 
meter measurements of the tides of the solid Earth (Will 1971d), and hence 
is no longer a viable theory. 


by four attributes: 

i) Curvature of Space-Geometry 

Metric theories of gravity predict that the spatial metrical properties 
of objects are describable by a non-Euclidean, or "curved" geometry, but 
differ in their predictions for the strength of this curvature. The amount 
of space - curvature which a standard mass produces in a given theory is mea- 
sured by a parameter 7 . 

ii) Nonlinearity of Gravity 

The superposition law for gravity in most theories is nonlinear (most 
theories predict that gravity itself also produces gravity), and the para- 
meter (3 measures the extent of nonlinearity of post-Newtonian gravity. Note 
that 7 and (3 are the usual Eddington-Robertson-Schif f parameters used to 
describe the "classical tests" of general relativity. 

iii) Preferred Universal Rest-Frame 

Some theories of gravity single out or produce a "preferred" frame, 
related to the mean rest-frame of the Universe. These theories have recently 
been shown to predict observable Solar-System effects related quantitatively 
to the Solar System's velocity through this frame [§ II, § III; see also 
Will (l971d), and Paper II in this series]. Three parameters Q^, CL and 
C suffice to measure the strength and specific manner in which motion in the 


2 





preferred frame affects the post-Newtonian metric. All Ct ' s are zero in 
theories which have no such "preferred-frame" effects (such as general 
relativity) . 

iv) Conservation Laws for Energy-Momentum 

Some metric gravitational theories do not possess a complete set of 
conserved energy-momentum quantities for isolated gravitating systems. A 
set of parameters £ measures the effects resulting from a breakdown of 
conservation laws. For point mass gravitational systems, two parameters 
and £ suffice in the post-Newtonian metric; for perfect- fluid systems, 
four are needed, , £ , and f^. More complicated models for matter may 

require additional ^'s to specify the complete post-Newtonian metric. Theories 
which do possess conservation laws (such as general relativity) predict that 
all £'s should be zero [see Will (1971c) for a discussion of conservation 
laws using the original PPN formalism]. 

Table 1 presents a summary of this new version of the PPN formalism. 

Each of the new parameters has an equivalent expression in the old forma- 
lisms of Nordtvedt (1968b) and Will (1971a), and these expressions are 
given in this table. 

Table 2 lists the values of these parameters for all viable metric 
theories which have been analysed to date (Ni 1972) and for one member of 
a new class of metric theories which are presented in this paper for the 
first time (§ III). 

In Table 3 the parameter dependences of a variety of experimental tests 
are presented using the new PPN formalism. Some of the effects listed in 
Table 3 are new and are discussed in Paper II of this series. 

This then, is the new PPN formalism in summary, and the remainder of 
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this paper is devoted to details. In § II we present and discuss the 
metric for the extended PPN formalism. This metric contains additional 
metric terms which are functions of the velocity of the observer's chosen 
coordinate system relative to some "preferred Universal rest- frame", and 
we show that these additional terms are necessary in order to handle 
"preferred- frame" theories of gravity from any frame of reference, more 
easily and consistently. Section III contains detailed discussions of 
two classes of "preferred- frame" theories: (a) theories which introduce 
the preferred frame as a basic postulate, called "stratified" theories 
(see Ni 1972), and (b) theories in which the preferred frame is established 
by cosmological fields (vector fields or second tensor fields) which are 
determined by the Universe as a whole. Concluding remarks are presented 
in § IV. An Appendix gives detailed calculation of the PPN metric of a 
vector-metric theory as an example of (b) above. 


II. THE EXTENDED PPN FORMALISM 

a) Key Features of the Old Formalism 
We first outline PPN formalism as described in Nordtvedt (1968b) and 

Will (1971a). The PPN framework has the following key features: 

/ 

i) Coordinate System : the formalism uses a coordinate system in which, 

at large distances from the matter, the metric becomes nearly Minkowskiian. 

ii) Matter Variables : gravity is produced by matter and energy, and 

the various matter variables which enter into the metric are: masses of 
particles (rru), rest-mass densities of fluid elements (p), specific internal 
energies (n), and pressures (p), as measured in local inertial frames comoving 
with the matter. The metric also contains positions and velocities (v) of 
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elements of matter as measured relative to the coordinate system. 


iii) PPN Parameters : 7 , p, Ct^, ^ 2 ’ £ 3 * ( see Table l). 


tach of the authors of this paper originally used his own set of PPN 
parameters. By using the parameters given in Table 1, we have united 
our separate notations, in addition to using parameters with more con- 
ceptual significance. Therefore, here we use the "new" parameters to 
review the old formalism. 

3 

iv) Metric of the Old PPN Formalism : 

3 

In this series of papers, Greek indices will take the values 1, 2, and 3; 
Roman indices will take the values 0, 1, 2, 3; and summation over repeated 
indices will be employed. Commas will denote partial differentiation, 
and semicolons will denote covariant differentiation. We will use units 
for which the velocity of light is unity and the Newtonian gravitational 
constant in the outer regions of the solar system, as measured today, is 
unity. 


A. 

g 00 


8 oa 


S Qp 


Perfect Fluid Metric. 

= 1 - 2U + 2£U 2 - (2y + 2 + a 3 + £ ) 3> 1 

- 2 [(3 7 - 2P + 1 + £ g ) 0 2 + (1 + £ 3 ) $ 3 + 3(7 + £ 4 )®jJ 

= '2(^-7 + 3 + - a 2 + £^) v a +-2(1+012- ^ , 

= - (i + 2 7 U) 5qj3 . 



(1) 
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Here U is the Newtonian gravitational potential given by 


U(x, t) = f 


p (x' , t) dx' 


the other potentials appearing in equations (l) are 




p* v' 2 dx' 


> ®2 = J 


p 1 U' dx' 


p' n' dx' 


. \ - ! 


p' dx' 


G = J 


p ' [ v 1 • (x - x')] 2 dx' 


p' v' dx' 

v a ~ J* ■ rx-TF T 


p* v* • (x - x' )(x - x' ) dx' 


The amount of rest-mass contained in a proper colume d?r of fluid is given 
by [see Will (1971a) for a discussion] 

pd^ = p(l + -| v 2 + 37U) dx 


B. Point-Mass Metric. 


“v / m k \ 

8 ~ =l - 2 l -J 


m k v t 


(27 + 1 + 0^ + £ x ) Z — — 

k k 


m m . 

2P + £ ) E — E 

k r k j/k r jk 


k- P 

+ ^1 J 3 ^~k ’ Hk^ ’ 

k r k 
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’occ 


5 op 


= |(^7 + 3 + ot x - a 2 + 5 X ) S ^ v k a + |( 1 + a 2 

k k 


m. 


(l + 2 7 Z S 

' k k' 


QP * 


?1> l “4 <X k • E k > 

k r k 


(5 con' t. ) 


a 


Here r, = x - x, , r . , = x . - x, and r, = r, , and 
~k ~ ~W ~ik ~x ~k k ~k ’ 


m, = T pd? 1 '" 

° k th body 


( 6 ) 


v) Choice of Gauge : an infinitesimal coordinate transformation 


/ i\t i i 
(x ) = x + gx , 


( 7 ) 


changes the metric to 


( 8 tJ ) - ®i j - 8*i. j - 6X J;1 


( 8 ) 


12 3 

We have chosen gx , §x , gx in order to make the spatial part of the metric 

(g ) diagonal and isotropic. We have chosen gx^ of the form 
Op 


0 

gx = A X Q , 


( 9 ) 


where A is an adjustable parameter, and where 


x(x, t) = - J p(x*,t) |x - x' | dx' [fluid] , 

= - Z m, |x - XjJ [point mass] , 

k 


( 10 ) 


in order that the metric component g ^ contain no term of the form 
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p' p" (x - x' ) • (x' - x" ) 

dx' dx" [fluid] , 

|x - x* | jx* - x" ~ ~ 

( 11 ) 

“k ~k ’ d ~k/ dt 
= - S — [point mass] 

k r k 

b) Lorentz Transformations and Preferred Frames 
Recent theoretical evidence has led us to conclude that the PPN metric 
in the form of equations (l) and (5) is valid for any metric theory of gravity 
only in a coordinate system which is at rest relative to a "preferred Univer- 
sal rest- frame". In fact, we will show (§ II. c) that there is a more general 
PPN metric which is valid in any chosen coordinate frame which moves at 
velocity w relative to the "Universal rest- frame ", and that the old PPN 
metric is a special case of this general metric corresponding to w = [ speed 
of chosen frame relative to Universe] = 0. Three pieces of evidence support 
this conclusion: 

i) A Lorentz transformation of the old PPN metric from one coordinate 
frame to another coordinate frame whose outer region moves with velocity 
w (|w|assumed small) relative to the first does not always yield the PPN 
metric in its initial form but introduces new metric functions which depend 
on w. Only when the PPN parameters satisfy the conditions 

0=1 = °, 0= 2 = 0 , 0= 3 = 0 , ( 12 ) 

do no such w- terms arise (Will 1971c, Nordtvedt 1969). (More precisely, 
when = °= 2 = 0^ = 0, any w- terms which do arise have no physical meaning; 
they can always be removed by an infinitesimal gauge transformation. See 
§ II. c for details.) 


= n 
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ii) A calculation by one of us (Will 1971d) showed that the Newtonian 
gravitational constant G as measured by a Cavendish experiment on Earth 
should depend on the Earth's velocity relative to some "preferred frame", 
for any theory whose PPN parameters violated equations (12). 

iii) A compilation and study by Ni (1972) of contemporary metric theories 
of gravity gives several specific examples of theories which single out such 
a "preferred Universal rest-frame". These theories all predict that at least 
one of the CK- parameters should be non-zero. We will discuss these theories 
in more detail in § III. a. 

There is additional motivation for using a more general PPN metric 
which is valid in any coordinate frame. Many computations of solar- system 
effects are simplified by working in the rest frame of the system being 
considered. Also, comparison of experimental results with the predictions 
of theory becomes more straightforward when the experimental data and the 
theoretical predictions are analysed in the same coordinate system. There - 
fore, for use in any frame moving at speed w relative to the Universal 
rest-frame, we will adopt a more general PPN metric which contains additional 
metric terms depending on w. Just what these additional terms are is spelled 
out in the next sub- section. 

c) The General PPN Metric 

We take as our starting point the PPN metric valid in the mean rest- 
frame of the Universe [eqs. (l) or (5)]. To obtain the form of the metric valid 
in a frame moving with velocity w relative to this frame, we apply a Lorentz 
transformation to the metric, assuming throughout that w is small, so that 
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The transformation from the 


our post-Newtonian expansion is still valid. 


4 


The post-Newtonian metric is an "expansion” in terms of small dimention- 
less quantities which characterize the gravitating system: Newtonian gravi- 

tational potential U ~ m^/r^, velocity v , specific internal energy II, 

— 6 

pressure relative to density p/p, all of which are < 10 in the solar system. 
In order to simplify the "bookkeeping" of the terms in the metric, one assigns 

them "orders of smallness" (Chandrasekhar 1965). In our notation U ~ m ^/ r j c 

2 

is 0(2), v is 0(l), v is 0(2) and so on. We will also assume that w is 0(l). 


Universal rest-frame coordinates (x, t) to the moving coordinates (l,r) has 
the form (to post-Newtonian order): 

. / i 12 3 4\ /_ 1 2\ i 

r=t(l+^w + g- w ) - x • w(l + -g w ) + 0(5) , 


12 1 

£ = X - (l + 7 T w ) wt + — (x • w) w + 0(4) 

~ ~ t— ~ cL ~ ~ ~ 


(13) 


The matter velocity in the two frames, v = dx/dt and v = d|/dx are related 
by 

v = v - w + 0(3) . (14) 

In order to express the metric functions in terms of (§, t) and v we make use 
of the expression given by Chandrasekhar and Contopoulous (1967) 


IS - V 


x - x' 


1 + ■ (s' 


w 


) (n* • v) + 0(4) 


(15) 


where 


n' = (x - x' )/|x - x' | 
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We use the standard transformation law (x° = t, = t). 


g, 4 (6 ,t) 


ij 


dx chc 

ait 


as 




( 16 ) 


along with equations (l), (4), (13), (l4), and (15) to obtain a metric which 
has the same functional form (in terms of t, v, etc.) as the old metric 
(in terms of x, t, v, etc.) except for additional terms which depend on w 
[we quote only the perfect- fluid case; an equivalent result is obtained using 
the point-mass metric eq. ( 5 )]: 


5 g oo = («! - a 2 - a 3 ) W u - (2a 3 - + a 2 + ^ - 1) w 

"ap ’ 


a 


CL 


+ (cl 2 + ” l) + OU w a \P U„ 


Q! 


68 0a 2 ^ a i " a 2 + ^1 ' w U + 2^2 ‘ ^1 + ^ w U oe » 


(17) 


88 ap 0 




where the potential U 


Cp 


U Qp X ,Qp + 8 OP 


is given by [cf. eq. ( 10 ) ] 

U , 


p(I',t)(S - i')„ (l - ?')„ 

-X— ^ , [fluid] 


£ - £ 


. |3 


= 2 


CL P 

m r r. 
k k k 


[ point mass ] 


(18) 


The introduction of the new variable w in the PPN metric also introduces 
an additional gauge freedom. By making a coordinate transformation of the 
form 
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( 19 ) 


t 1 / „ , \ cc 

T T 2 °2 + ^ 1 “ 1 ) w > 

! + = ! > 

and using equation (8), we can transform the metric into a gauge in which 

Ct 

no term of the form w W^, appears. Note that this gauge transformation is 
completely independent of the gauge change [eq. (9)] which led to the form 
of the old PPN metric in equation (l). The final result for the new PPN 
metric for perfect fluids and for point masses, is given in Table It-. 

Table 4 shows clearly the usefulness of the new parameters CC^, QL^, and 
QL. These parameters label the terms in the metric which depend on w and 
indicate, for a given theory, how large and what kind of solar-system effects 
will be produced by motion relative to the Universal rest-frame . Paper II 
in this series (Nordtvedt and Will 1972) will catalogue such solar- system 
effects and present evidence to put experimental limits on the a' s, thereby 
ruling out many theories of gravity which were thought at one time to be 
viable. 

We close this section with a recipe for computing the PPN parameter 
values for any metric theory of gravity: 

i) Pick a coordinate system whose metric becomes Minkowskiian far 
from the local distribution of matter. It doesn't matter whether the coor- 
dinate system is moving or at rest relative to the Universe, although for 
many theories (cf. § III, and Appendix) calculations are simplified using 
the "rest" coordinate frame of the Universe. 

ii) Use the standard techniques of Chandrasekhar (1965) or Einstein, 
Infeld, and Hoffman (1938) to compute the post-Newtonian metric of the matter 
in this coordinate system. 
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iii) Compare this metric term by term with the equations in Table 4 
and read off the values of the PPN parameters. Even if you used the "rest" 
coordinate frame, in which w = 0, you will still be able to obtain values for 
all the PPN parameters. 

iv) To obtain the metric valid in any chosen frame, use the equations 
in Table 4, along with the PPN parameter values just computed, the velocity 
w of the chosen frame relative to the Universe, and the velocities v of 
matter relative to the chosen frame (the other variables: mass, density, 
etc., are all determined in frames comoving with each element of matter and 
are unaffected by this change in coordinate frame). 

III. THEORIES OF GRAVITY WITH A PREFERRED 
UNIVERSAL RE ST- FRAME 

There is a variety of gravitational theories which single out some 
"preferred frame" in the Universe, and which have been shown to predict 
observable effects due to motion through this frame. These theories fall 
into two classes: (a) Theories which postulate the existence of a preferred 
frame directly and incorporate it into their calculational rules. These 
theories endow the Universe with "prior geometry", and use this prior- 
geometric structure as a foundation for the theory.^ In § III. a we discuss 

5 

For a detailed discussion of "prior- geometry" see Misner, Thorne, and 
Wheeler (1972). 


a particularly simple subclass of these theories, known as "stratified 

i 

theories with time -orthogonal conformally flat space slices" (Ni 1972). 
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(b) Theories which assume no "prior geometry", but which possess, in addition 
to the metric, vector or tensor cosmological fields whose source is the matter 
in the Universe. In these theories, the particular structure of the vector 
and tensor fields determines (in a Machian sense) the "preferred frame". 

These theories are examined in § Ill.b and in an Appendix. 

a) Stratified Theories With Time-Orthogonal 
Conformally Flat Space Slices 

This is a particulary simple subclass of theories of gravity which 

postulate a "prior geometry". These theories are devised using the following 

0 

prescription: The Universe's large-scale distribution of matter determines 

0 

This discussion is adapted from Ni (1972). 


a preferred reference frame whose space slices ("strata") are conformally 
flat, although the full spacetime is not. In this preferred frame, the 
metric has the form 

ds 2 = e 2f ^ dt 2 - e 2g ^ (dx 2 + dy 2 + dz 2 ) , (2o) 


where cp is a scalar field. In geometric, coordinate-free language, such 
theories have (i) a background, flat metric (ii) a Universal time coor- 

dinate t (a scalar field) which is convariantly constant and has timelike 
gradients with respect to (iii) a scalar gravitational field cp; and (iv) 
a metric g constructed from t, and cp by 

g = e 28(f)^ + [e 2*(<p) . e 2g(-f)]dt g dt . 


( 21 ) 


These theories differ from one another by their choice of the function f(cp) 



and g(cp) and by their field equations for cp. 

Ni (1972) has analysed theories of this type, whose authors include: 
Papapetrou (1954a, b,c), Yilmaz (1958, 1962), Page and Tupper (1968), and 
Rosen (1971a, b). Two additional theories were devised by Ni himself: a 
Lagrangian Stratified Theory and a General Stratified Theory (Ni 1972). A 
special case of Ni's "General Stratified Theory" is Coleman's (1971) theory. 
Einstein (1912) ( not general relativity.') and Whitrow and Morduch (i960, 

1965) devised stratified theories of gravity, but their theories disagree 
violently with light deflection and time-delay experiments, and so will not 
be considered further. 

We will not discuss any of these theories individually; the reader is 
referred to Ni (1972) for details. We will only present Ni's results, in 
the form of PPN parameter values for each theory (Table 2). For complete- 
ness, Table 2 also gives the PPN parameter values for general relativity 
and for scalar-tensor theories [generalizations of Brans-Dicke- Jordan theory; 
see Ni (1972) for details and references]. We conclude this subsection with 
a theorem concerning the PPN parameter o:^ for stratified theories: 

In every stratified theory of gravity with time-orthogonal, conformally 
flat space slices, the PPN parameter has the value 

a ± = - 4(1 + y) . ( 22 ) 

In order to agree with light deflection and radar time delay experiments 
[7 ~ 1; see Thorne, Will, and Ni (1971) for a review], these theories must 
therefore have 

~ - 8 . (23) 
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Note that all the stratified theories in Table 2 obey equation (24). In 
Paper II of this series, we will make use of this result to show that all 
the stratified theories in Table 2 disagree violently with experiment . 

The proof of this theorem goes as follows: Pick any stratified theory 
with time-orthogonal conformally flat space slices. Using a coordinate 
system at rest with respect to the Universe, compute the post- Newtonian 
metric due to an arbitrary configuration of matter. To put this metric 
into the "standard" PPN gauge, it may be necessary to apply an infinitesimal 
gauge transformation [eq. (9)] 


t + = t + A x Q , x + = x , (24) 

for some value of A. Since g^ was initially identically zero [by assump- 
tion, cf. eq. (20)], in the new gauge it becomes 


( g ocP A x , oa 


= - A v a + A w a 


(25) 


By comparing this with the PPN metric in Table 4 we obtain 

-^(by + 3+a 1 -a 2 + £^) = -A , 
i(l *<* 2 - Cj) - * • 

Adding equations (26) and (27) we obtain, finally, 

Oj_ = - (bj + 4) , 

independent of gauge, A. Q.E.D. 


(26) 

(27) 


( 28 ) 
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b) Lagrangian-Based Metric Theories with Additional 
Vector and Tensor Fields 


Here we take an approach which is different then that of § III. a, and 
analyse preferred- frame theories of gravity based on a Lagrangian formula- 
tion which initially assumes no preferred frames (no "prior geometry"). 

Since gravitation appears to be a long-range, purely attractive force, 
one might expect the Universe's global matter distribution to affect local 
gravitational physics, and to play a dominant role in establishing any 
"preferred" reference frame. In fact, from this "Machian" point of view, 
the mystery seems to be: How can a metric theory of gravity avoid at all 
having a preferred frame related to the Universe rest- frame ? General 
relativity and the Brans-Dicke theory obviously avoid preferred frame 
effects (cf. their parameter values = 0 in Table 2). Before 

we create theories which dx> have preferred frames, we must first shed some 
light on how these exceptional theories avoid them. 

Consider a local gravitational system, such as the solar system, which 
is embedded in the Universe. We separate the computation of the metric into 
two parts: a Universal or cosmological solution, and a "local" solution. 

From this viewpoint, the Universe affects the local gravitational physics of 
the system by establishing the boundary conditions (at a boundary "far" from 
the matter) for the various fields generated by the local system. The local 
system "feels" its relationship to the Universe via the asymptotic field 


values of the fields present; metric g^, scalar field cp, vector field K^, 
second-rank tensor field C^, and so on. Several conclusions follow: 

l) A theory which contains solely a metric field yields local gravi- 


tational physics which is identical in all asymptotic Lorentz frames, and 


17 








which does not change with Universe evolution. In addition, the Newtonian 


gravitational constant G is unaffected by the proximity of matter . All this 
follows from the invariance properties of t)^ (the asymptotic form of g^), 
the only field coupling the local system asymptotically to the Universe, 
and from general covariance, which allows us to find a coordinate system in 
which the metric field takes this Minkowskii form at the boundary between 
the Universe and the local system. 

2) A theory which contains a metric field and a scalar field cp yields 


)hysics which is identical in all asymptotic Lorentz frames, but which may 


vary with Universe evolution; G may be affected by proximity of matter. 


These conclusions follow from invariance of both -n . . and cp under Lorentz 

'ij T 

transformations, but now cp may vary with Universe evolution and may depend 
on the proximity of matter. 

3) A theory which contains a vector field K. and/or an additional 


second- rank tensor field C. . yields local physics which does depend on 


— J 

motion relative to a preferred Universe rest- frame, and which may vary 


with evolution of the Universe: G may be affected by proximity of matter. 


This follows because the asymptotic values of K^, and are not invariant 
under Lorentz transformations (an exception would be [ aa yiUF LwL ^ j propor- 

tional to n . . ) . 

ij 

In summary, it is the Lorentz invariance of the asymptotic fields cp 
and which makes it impossible to have preferred- frame effects in theories 
containing solely those fields. Thus we must appeal to vector- or tensor- 
metric theories for preferred- frame effects. 

It is commonly believed that cosmological vector fields and additional 
second-rank tensor fields are absent from physics. The Hughes-Drever 
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experiment (measurement of the isotropy of inertial mass) and a variety of 
laboratory "ether-drift" experiments have been used by Dicke [ (1964); see 
also Peebles and Dicke (1962); Peebles (1962)] to rule out these fields. 
Closer examination of these arguments shows that the experimental evidence 
rules out only those vector and "second tensor" cosmological fields which 
couple directly to matter . Since these experiments were performed under 
conditions where the effects of gravity were negligible, they do not rule 
out vector and second- tensor fields which couple only to gravity . 

We now proceed to devise such theories. We assume these theories are 
derivable from a coordinate- invariant Lagrangian, with no a priori assump- 
tion of a special coordinate frame or of any prescribed form for the tensors 
in the theory. All physical quantities (metric, vector fields, tensor fields, 
matter variables) must be calculated in terms of other physical quantities. 

The Lagrangian for such theories has the general form 


L = L m (q,p,g..) 


+ L ( g . g . . . # K. , K. . , C . C , •••) , 

g ij 3-J, k' i’ 1, y xj’ i J, k' 


(29) 


where q and p symbolically represent the dynamical variables of matter, g 


ij 


is the metric and K., C. 

i ij 


. are possible cosmological vector and tensor 
fields. The matter Lagrangian is related to the stress-energy tensor for 
matter and non- gravitational fields by the functional derivative 


1 

/-g 



(As I J 



( 30 ) 


The additional fields in the Lagrangian do not couple directly to matter; 
this leads to the equation of motion for matter 
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= 0 


(31) 



> 


which yields geodesic motion for test bodies. This form for the Lagrangian 
thus guarantees that the theory satisfy the postulates of a "metric" theory 
[for a discussion of these postulates see Thorne and Will (1971)]. 

For simplicity, we restrict ourselves to theories which have linear 
field equations for all the supplementary fields; nevertheless, a wide 
variety of theories is possible. For vector-metric theories, the general 
Lagrangian has the form 




V. - L „ + k l R + k 2 K 1 K R + "s K i K J R 

+ ^ K 1;j + k 5 K K J;l + \ K ,' 1 K.^J , 


(32) 


where R. . and R are the Ricci tensor and scalar formed from the metric, k, , 
ij ’ r 

. .., kg are dimensionless coupling constants. Tensor-metric theories have 


even more structure (we assume is symmetric): 


L m “ L m + V R + V c u clJ R + V< C ij ^ 




R 


+ k,, • C . . 
4 U 


= lj rW - V ^ RiJ 


+ k ’ C. . C lj;k + k’ C. . C lk;j + k ' C.. V j 

6 ij;k 7 ij;k 8 ik; 


n 


(33) 


+ k ' jk£ ;i , k ij Im 

9 C ii; 8 V + k 10 C ij;k C im; 8 8 


We point out here that the above equation (33) alone does not reveal the 
full richness of possible Lagrangian terms for tensor-metric theories. If 
is a non-singular field, one can define Christoffel symbols in terms 
of C . Also, third-rank tensors S 1 can be produced by taking the difference 
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of the Christoffel symbols 


jk 


= r 


jk 


(g) - 


jk 


CC) 


Covariant derivatives can be defined with respect to each tensor, or 

and tensor densities can be formed using either 

or 


/|det C | * 


All these quantities can be used to construct Lagrangian terms. 

A future paper in this series will give a complete study of these theories. 
Such a study will discuss consistency of the fields equations, the post- 
Newtonian limit and the PPN parameter values, cosmological solutions, and 
positive-definiteness of free-field energies. 

To show that theories of this type can produce apparent agreement with 
the "three classical tests" of general relativity, yet also possess preferred- 
frame effects, we now quote the PPN parameter values for a simple illustra- 
tive case of a vector-metric theory. Detailed calculations have been rele- 
gated to an Appendix. This theory has as its action integral 

A = J d**x /- g [G 0 t m + E + K.. . g lk gJ'] (Zb) 

where is the unrenormalized gravitational constant. This theory corres- 
ponds to the case [eq. (32)] = k^ = G Q = k^ = = k g = 0. 

The calculations in the Appendix yield a PPN metric with parameter 
values (see also Table 2) 
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7 = 0=1 


°1 = a 3 = h = ^2 = ^3 = ^ = 0 (35) 

a 2 = £/( 1 + \ K 2 ) . 

2 

Here K is the squared magnitude of the vector field, evaluated far from 
the matter, i.e., 

K 2 = [g lj K K ] . (36) 

- 1 far from matter 

The gravitational constant is renormalized in this theory, and so changes 
in time as the Universe evolves: 




(37) 


IV. CONCLUSIONS 

In this paper, we have presented a new version of the Parametrized 
Post-Newtonian Formalism, and have discussed preferred- frame theories of 
gravity. In Paper II of this series we will use this new formalism to com- 
pute observable solar-system effects predicted by preferred- frame theories 
and will present experimental evidence to rule many of these theories out. 

The authors wish to thank the National Science Foundation and the 
Montana State University Physics Department for their support and hospitality 
during the 1971 N. S.F. Summer Workshop in Selected Topics in Theoretical 
Physics (GZ-1919-NSF), where the initial research for this paper was carried 
out. 

We also thank Ron Hellings, Wei-Tou Ni, and Kip S. Thorne for helpful 
discussions. 
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APPENDIX 


A VECTOR-METRIC THEORY OF GRAVITY 

Here we discuss in detail the vector-metric theory described in § Ill.b. 
This theory has as its action integral [cf. eq. (34)] 


a= 8 [g 0 L m+ R + K 1;J g lk gj'] , 

where L is the matter Lagrangian, R is the scalar curvature formed from 


(Al) 


the metric g^y is the cosmological vector field and Gq is the (unrenor- 
malized) gravitational constant. By varying the action with respect to g^ 
and in the usual way, we obtain the field equations 


K. 


J 


= 0 


R. . — "o R § . . — 8 it G T. . + 0. . 
xj 2 6 xj 0 xj ij 


(A2) 

(A3) 


where R. . is the Ricci tensor, T. . is the stress-energy tensor for matter 
ij ' xj 

and non-gravitational fields, and 0. . is the vector-field stress-energy 

^ J 

tensor given by 


0 . . 

xj 


K/ m + K . K 111 . - 3 K™ K ;n g.. 
x;m j m;x ;j 2 ;n m °xj 

S. . - K. S m . - K. S m ) 

2 \ xj x j J x';m 


(A4) 


where 


S. = K. . + K. . 
xj x;j j;x 


(A5) 


It is straightforward, though tedious, to verify using equations (A2), 
(A4), and (A5) that 


0 * J = 0 
xj ’ 


(A6) 


and thus that 
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(A7) 


We now use these equations to compute the metric produced by a con- 
figuration of perfect fluid, in the post-Newtonian approximation, using the 
standard techniques of Chandrasekhar (1965). [Equivalent results can be 
obtained for a system of point masses using the standard EIH (Einstein, 
Infeld, and Hoffman 1938) method. ] The metric is expanded about the flat- 
space Minkowskii metric according to 


g..(x, t) = Ti. . + h..(x, t) 
& ij ~ 'ij ij 


(A8) 


where, far from the matter, h^ tends asymptotically toward zero. We also 
expand the vector field about its asymptotic value. However, in order 
to simplify the calculation, we will work in a coordinate system in which 
the asymptotic vector field has only a time component; i.e. asymptotically 


K a = 0 . (A9) 

Then, in this coordinate system, we have 

|K| 2 5 q ij K. K. = K q 2 , (A10) 

far from the configuration of fluid. The frame in which = K and K = 0 

& 0 a 

is presumably the rest-frame of the Universe's smoothed_out distribution of 
matter. The expansion of the vector field can thus be written 


Kq = K + <p(x, t) , (All) 

K a =k a ( *,t) , (A12) 

where cp is 0(2) and k is 0(3) (cf. footnote k), and both go to zero far 


2 k 



from the fluid. 


The perfect-fluid stress-energy tensor is given by (cf. Will 1971a) 


T ij = (P + P n + p) u i u j ~ P g t j > 


(A13) 


where u^ is the four- velocity of the fluid. The field equation (A3) can 
be written in the form 


R u = 8 ” G o (T u - i 1 hP + ( e ij - 1 e *u> ’ 


(Alh- ) 


where T and 0 are the traces of T. . and 0. . respectively. 

ij , ij 

We wish to obtain h^ 0 to O(b-), h^ to 0(3), and h to 0(2). The 
calculation proceeds along the following lines: 

i) Calculate h^ to 0(2). This corresponds to the theory's Newtonian 
limit. To the required order, we have 


^00 2 ^ ^ nn ’ ^ nn ^ ^ 


0 


00 


Ak 2 - 2 


00 ’ 00 
v h oo ’ 


e = - k r 


00 


(A15) 


Equation (Alb-) then becomes 


V 2 h Q0 = 8 * pjc^d + -| K 2 )] , 


whose solution is 


h oo 


2u[g 0 /(1+A k 2 )] . 


(A16) 


(A17) 


Because we have chosen units so that the gravitational constant as measured 
far from the solar system, today is unity, we have 

G, 


0 


1 + 


4 k 2 


= G = 1 [ today ] 


(A18) 
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Note that if K varies slowly as the Universe evolves, then the gravitational 
constant G may vary from its present value with time at a rate 


l/dG\ K /dK\ 

G\dt/ " 1 v 2\dt) * 

' 1 + K ' ' 


(A19) 


ii) Calculate h^ to 0(2). By making a particular choice of gauge: 


h i,op " h7 p,07 “ h7 a,P7 ° ' 


(A20) 


we give the field equation for (eq. AlU) the simple form (to the required 

order) 


whose solution is 


V h QJ3 = 8 * P S OP ' 


h ap 2 u 5 ap 


(A21) 


(A22) 


iii) Calculate cp and k to 0(2) and 0(3). The field equation (A2) for 
i = 0 becomes, to the necessary order 


V (cp + K U) = 0 


(A23) 


Since cp must vanish far from the fluid, we have 


cp = - K U 


(A24) 


If we choose the gauge in which 


1 uP 

2 


hP p,o- hP o,p=' u ,ol K 2 /( 1 + l K2 )] - 


(A25) 


Then equation (A2) for i = OL becomes 


v 2 


k “^ Kh “ + H + ^b < » 


= 0 , 


(A26) 
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where X is given by equation (10). Equation (A26) has the solution 


k a4 Kh oa-¥ K ( 1+ 771^) x ,oa * 


(A27) 


iv) Calculate h^ to 0(3). With the choice of gauge (eq. A25), the 
field equation (AlU) simplifies to 


y2 h oa = ' 16 p v a - 1(1 - 


K 2 


X 


or 


1 + I K 2 / '°“ 


(A28) 




a ' 2 V"l + |K 2 J % ’ <a ' 


(A29) 


Using the fact that 


X ,ca V oc ~ w a * 


(A30) 


we get 


K 2 


OCC 2 


1 - + A 1 * 


K 2 


1 + \ K 2 , 


W. 


a 


(A31) 


v) Calculate h^ to O(^). To the appropriate order, equation (AlU) 


along with equations (A13), (A17), (A2^), (A25), and (A27) gives 


V 2 


00 


8 it p + 2 V 2 (U 2 ) + 16jtpv 2 + 8rtpn 


(A32) 


0)l 


PZi tT 


whose solution is 

(A33) 

The final form of the metric is 


00 


= - 2 U + 2 U - It - 4 &, 


- 2 *3 - 


6 $, 
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? 0Q = l- 2U+2U 2 -i4-<2> 1 -ll-$ 2 -2$ 3 -6$^ 


g oa 2 V " 


A? H + K 1+ 77 


1 k 2 / W a 


(AS4) 


g Qp - (1 + 2 u) 5 Qp 


By comparing equation (A3h) with the metric in Table k, we can read off the 
values of the PPN parameters, as listed in Table 2. 

Thus this theory does predict that local gravitational physics should 
depend on motion relative to the preferred Universal rest- frame (since 
0), and should vary with evolution of the Universe [since K may vary with 
time; cf. eq. (A19)]. Notice that, in a Robertson-Walker universe with 
expansion parameter a(t), evolves according to ("dot" denotes d/dt) 


K q + 3 K Q (a/a) - 3(a/a ) 2 K Q = 0 


K = 0 

a 


(A35) 


However, this theory predicts that the local gravitational constant G 
should be unaffected by nearby matter. In the presence of matter (and in 
the Universal rest- frame), this locally measured G is given by [see Nordtvedt 
(1970b), Will ( 1971d) ] 

G local = 1 - (4 P - r - 3 - ?2 ) U external . (A36) 

Since (3 = 7 = 1, and = 0 for this theory, we find G^ oca ^ *- s unaffected by 
nearby matter. This can be seen in another way, as follows: Consider a 

gravitating system in the presence of an external mass, whose Newtonian 
gravitational potential is U . When the local system and the external mass 

cXt 
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are analysed together in an asymptotically flat coordinate system, the 
metric and the vector field are given by, to first order [ cf . eqs. (A2U), 
(A34) ] 


g 00 = 1 - 2 u - 2 U e 

K 0 = K < 1 - u - U eX t> 


(A37) 


The gravitational constant G is determined by the value of K Q far from the 
system; i.e. [ cf . eq. (A17)] 



(A38) 


We now look at the local gravitational system itself using a coordinate 
system which becomes asymptotically flat (to some desired precision) far 
from the local system but between the local system and the external mass. 
The transformation from the first coordinate system to this one has the 
form 


which gives 


t 

new 


t old^ 1 ' 


U ) 
ext' 


^ 8 0(Pnew 1 " 2 u 


(A3 9) 


(Vnew 


K( 1 - U) 


(AkO) 


Thus the transformation which removes the external field from the metric also 
removes the external field from the vector field. Therefore the gravitational 
constant as measured in the local system is unaffected by the external mass. 
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THE NEW PPN PARAMETERS AND THEIR INTERPRETATION 
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Adapted from Thome, Will, and Ni (1971). 

The metric in the old formalism was valid only in a coordinate frame at rest in the Universe 
rest- frame (see §Il). 



METRIC THEORIES OF GRAVITY AND THEIR PPN PARAMETER VALUES 
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Ni (1972) for detailed discussion of all the theories in this table except the Vector-Metric Theory 
: § Ill.b, Appendix), and Coleman's (1971) theory. 
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TABLE 3 


PPN PARAMETER DEPENDENCE OF VARIOUS EXPERIMENTAL TESTS 


Experimental Test or 
Observable Effect 


PPN Parameter 
Dependence 


Light Deflection, Radar 
Time Delay 

Perihelion Shift 

a) Classical 

. * 

b) Due to Motion through 

Universal rest- frame 
Geodetic Precession of a gyroscope 
Dragging of Inertial Frames 

Equivalence Principle Breakdown 

(m . 4 m. . . , ) 

passive ' inertial 

a) Isotropic 

b) Anisotropic 

Perturbations in Earthbound 
Gravimeter Measurements: 
a) Variation in G 

i) Due to field of Sun 
and Planets 

ii) Due to motion through^ 
Universal rest- frame 
iii) Anisotropy due to motion 
through rest- frame 


■|(i + y) 

|(2 + 2y - p) 

y ^2* 

|(i + 27) 

|(1 + 7 + { 0^) 


- (Up - 7 _ 3 - 0^ + a 2 - 5j_) 

- (a 2 + ^2 - Si) 


Up - 7 - 3 - 
a 2 + °3 - a i 
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TABLE 3 (continued) 


Experimental Test or 
Observable Effect 


PPN Parameter 
Dependence 

b) Other Effects 

i) Due to external field 
gradients 

“l 


ii) Due to internal structure 
of the Earth 

bp - 7 

- 3 - a 2 + $i - 2^ 2 

Periodic Perturbations in orbits 

due to Motion through Universe 

* 

re st- frame 

®1> ^2? 

a 3 


See Nordtvedt and Will (1972) (Paper II). 

^See Will (l971d), Nordtvedt and Will (1972) (Paper II). 
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TABLE 4 


THE PPN METRIC VALID IN ANY FRAME 


Perfect-Fluid Metric 

g 00 = 1 - 2U + 2PU 2 - (27 + 2 + a 

- 2 [(3/ - 2f3 + 1 + £ 2 ) $ 2 + (1 + ? 3 ) $ 3 + 3(7 + ®jJ 
+ (a^ - a 2 - a^) w 2 u + a 2 w a - ( 2 a 3 - a^) 

s 0a = 2^ 7 + 3 + a i ~ a 2 + ?i) v a + H^ 1 + a 2 ' Sf) W a 
+ -|(o: 1 - 2 a 2 ) vPSj + a 2 w P , 

g ap = ■ ( 1 + 2 ? u ) 5 ap 


Point-Mass Metric 


“k 


/ “vX “v m - 

= 1-2 E -= + 2p E - 2(1 - 2p + E E -J- 

k r k v k k' 2 k k #k jk 

m v m, p 

- (27 + 1 + a + £ — — + $! s — (v k • r k ) 

k k k r, 

k 

m, m, 2 tn, w • v, 

+ (a - a - a ) W 2 e ~r + oc £ - 5 f- (w • O - ( 2 «, - a,) £ 7 — — 

123 k r k 2 k r, k 3i k r k 


’oa 2 


i(47 + 3 + a x - a 2 + s 


a 

m v. - m 

k * + |(X + a 2 - ?1 ) 5 -^(v k -r k )r k “ 


k r k 


k r. 


m, 

+ 2 ( Ql - 2a,,) » E — + a 2 E 7% (W • r fc ) r k 

k k 


a 


k r. 


"(> 


5 0£ l 1 + 27 I rj 5 ap 
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I. INTRODUCTION AND SUMMARY 


In Paper I of this series (Will and Nordtvedt 1972) we presented a 
new version of the Parametrized Post- Newtonian (PPN) Formalism, the theo- 
retical tool used to study experimental tests of metric gravitational theories. 
An important feature of this new version of the PPN Formalism is a simple 
and consistent prescription for analysing "preferred-frame" theories of 
gravity. By "preferred- frame" theories we mean metric theories which single 
out the frame which comoves with the Universe's smoothed-out distribution of 
matter as a "preferred" reference frame. Several such theories are known; 
those studied to date fall into two classes: i) "Stratified Theories" (with 
"time-orthogonal conformally flat space-slices") whose authors include Page 
and Tupper, Yilmaz, Papapetrou, Rosen, Ni, and Coleman; and ii) theories with 
additional cosmological vector or tensor fields; one of these, a vector-metric 
theory was analysed in Paper I. 

Even though these theories single out a preferred Universal rest-frame, 
they have until recently been considered to be viable alternatives to general 
relativity, since their predictions for the "standard tests" — gravitational 
redshift, light bending, radar time delay, and perihelion shift — are all in 
agreement with observations. But a closer examination of these theoretical 
predictions reveals that they agree with observations only if the solar system 
is assumed to be at rest relative to the preferred Universal rest-frame . 

A much more reasonable assumption is that the solar system moves through 
the Universe with a velocity of around 200 km/sec, resulting from its orbital 
motion around the galaxy and from the peculiar motion of the galaxy itself. 
Measurements of anisotropies in galactic redshifts and in the cosmic micro- 
wave radiation lend some support to this assumption (see Sciama [1971] for 


/ 
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a review of the relevant data; see also §IV). In this paper we will show 
that such a 200 km/sec motion produces, according to preferred- frame theories 
of gravity, observable solar-system effects which are in disagreement with 
observation. 

Throughout this paper, we will use the language and notation of the 
new PPN formalism discussed in Paper I of this series; the reader is urged 
to consult this paper for details. This formalism treats the post-Newtonian 
limit of metric theories of gravity in terms of a series of metric parameters, 
whose values vary from theory to theory. These PPN parameters have a physical 
or conceptual significance; in particular, the three parameters OC CC^, and 
measure to what extent a theory of gravity singles out a "preferred uni- 
versal rest- frame". This can be seen from the form of the PPN metric in 
Table k of Paper I. When written in a coordinate frame which moves with velo- 
city w relative to some "preferred frame", the PPN metric may contain terms 
which depend on w; these additional terms are (cf. Table k, Paper i)^ 


The units, notation and conventions used in this paper are the same as 
those used in Paper I. 


5S 00 = (“i - - Qj) u + a 2 W a J 3 Uqp - (2^ - a,) »“ v Q , 

S8 oct " ! (a l ' S3 2 ) u + a 2 " P u ae < 


58 ap 0 ' 


( 1 ) 


where U is the Newtonian gravitational potential. The potentials and 
are defined in Paper I. 

These additional metric terms may produce observable solar- system effects. 
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which depend on the PPN parameters Ot , QL , and QL, and on our velocity w 
relative to the Universe. In this paper we derive and catalogue these 
"preferred- frame" effects, and write them in the form 


( size of\ 
effect / 


\ 


/ 


/combination of 
a 2 , and a 


YV 


( amplitude, which\ 
depends on w J 


( 2 ) 


In this form, the predicted sizes of "preferred- frame" effects can be com- 
pared with experimental or observational data and limits put on the values 
of the parameters a^, CCg, and 0! , once we assume a value for w. In this 
paper, we assume that the solar system 1 s motion through the Universe is due 
to its (nearly circular) orbital motion around the Galaxy, i.e. we assume 
w is 200 km/sec in the direction given by galactic coordinates = 90°, 



The "preferred-frame" effects which we have studied fall into two 
classes: geophysical (effects which alter the structure of the Earth) and 
orbital (effects which alter the orbital motions of planets and satellites). 

a) Geophysical Effects 

One of the authors (Will 1971) has previously used the PPN formalism 
to show that the Newtonian gravitational constant G as measured by means 
of Cavendish experiments may depend on the observer's velocity relative to 
a Universal rest-frame. We have confirmed and generalized this result using 
the new version of the PPN formalism discussed in Paper I: by considering 
the Earth' s gravitational force on a gravimeter at rest on the surface of 
the Earth, we obtain for the measured value of G, 
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( 3 ) 



where M is the mass of the Earthy R its radius and I its spherical moment 
of inertia: w is the Earth's velocity through the preferred Universal rest- 
frame, and e^ is a unit vector joining the gravimeter and the center of the 
Earth. The Earth's velocity w is made up of two parts, a uniform velocity 
w of the solar system relative to the preferred frame, and the Earth's 
orbital velocity v around the Sun, thus 


2 2 _ 2 
w_ = w + 2w • v + v 
© ~ 



e ) 2 = (w 
~r 


e ) + 2(w • e ) (v • 


e ) 
r 


\2 
e ) 
~r 


00 


So because of the Earth's rotation (changing e^) and orbital motion (changing 

v), there will be variations in the gravimeter measurements of G, given by 

-9 

(we retain only terms which vary with amplitudes larger than 10 G) 


AG/G - (| 


a 2 + 





V 



e r ) 2 + 2(w 





( 5 ) 


where we have used the fact that, for the Earth, 

I - | MR 2 . (6) 

The most pronounced effect is the anisotropy in G produced by the 
(w • e ) term (eq.[3]). This anisotropy has been shown (Will 1971) to 

vI7 ~ a" 

produce "tides of the solid Earth", i.e. variations in the acceleration 
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g measured by a gravimeter, which are completely analogous to the tides 
produced by the Moon and Sun. The most important of these Earth- tides is 
a 12-hour sidereal- time tide, produced by the (w • e ) term of equation 
(5), with amplitude 

(A g / g ) ppN - « 2 (3 X 10" 8 ) cos 2 (L) , (7) 

where L is the latitude of the gravimeter's location on the Earth. Will 
(1971) has examined Earth- tide data and found that any discrepancy between 

Newtonian theory and experiment for this component of the tides must be 

9 

less than one part in 10 at any latitude for which reliable gravimeter 
data is available. By comparing equation (7) with this experimental limit, 
we find that Oi^ must satisfy 

\a 2 \ < 3 x 10 " 2 . (8) 

In §11 we summarize other "Earth- tides" of various frequencies produced 
by the anisotropic terms in equation (5). 

As the Earth orbit the Sun, the (w • v) and (w • e f ) (v • e^) terms 
in equation (5) vary with a period of one sidereal year. Because the 
Earth is a gravitationally bound object, this variation in G causes the 
Earth to expand and contract spherically. This "breathing" of the Earth 
results in a changing moment of inertia, which in turn (by conservation of 
angular momentum) causes a yearly variation in the Earth's rotation fre- 
quency Q with amplitude, according to the PFN formalism (cf. §Il) 

(tfl/fl) ppN « (| a 2 + a 3 - aj (3 x lo" 9 ) . (9) 

The observed yearly variation in the Earth's rotation rate (measured by 
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comparing astronomical time with atomic time standards) has an amplitude 
of 

(^fi/a) 0BSERVE D « ^ X 10 (10) 

(Smith and Tucker 1953). But this variation can be readily understood using 
Newtonian geophysics: it is produced by an annual variation in the angular 
momentum of the atmosphere due to seasonal changes in wind patterns and by 
a long-period (one year) Earth- tide produced by the Sun (Mintz and Munk 
1953; see also Melchior 1966). These calculations yield agreement with the 
observed variation in 0 with uncertainties around 15 percent, hence the PPN 
variation in the Earth's rotation rate must satisfy 

(#) ppN < 6 X 10' 10 . (11) 

Equations (9) and (ll) thus show that the PPN parameter combination 
(2/3 a 2 + a 3 - ) must satisfy 

| a 2 + a 3 - aj < 0.2 . (12) 

In §11 we discuss other, smaller variations in the Earth's rotation 
rate produced by the velocity-dependent terms in equation (5). 

b) Orbital Effects 

The solar system' s 200 km/sec motion through the Universe may produce, 
according to the PPN formalism a variety of observable effects, both secular 
and periodic, in the orbits of planets. The most important of these effects 
is an anomalous perihelion shift for the planets, given by (in radians 
per orbit) 
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( 13 ) 


AjJ 


2 „ = | (ar + 2 - p) 


3M 


1 r/ /M\ l/2 W Q 1 / 2 

" ¥ l\pj e ~ 8 a 2 W P 


w n 2 ) 


+ A a JM 

+ 2 a 3 M 




where M is the mass of the Sun; p and e 0) are the semi-latus rectum 
and eccentricity of the planet's orbit; Wp and w^ are the components of 
w in the plane of the planet's orbit, Wp in the direction of the perihelion, 
Wq at right angles to Wp; |ft| is the gravitational self-energy of the Sun 
and A is the Sun's rotational angular velocity. The first term in equation 
(13) is the "classical" perihelion shift, which depends on the PPN para- 
meters 7 and 8 (see Paper i) and which would be present even if the solar 
system were at rest in the Universe (w = 0). For Mercury (?) and the 
Earth (©), equation (13) yields, in seconds of arc per century (§IIl): 


(43^ > PP N = + a - P)] + 35 Oj + 8 

= ‘‘l^ 27 7 2 - ?)| 7 57 “l 


0C 2 - 4 X 10 0^ , 


+ a 2 - 7 x 10 a 3 


(14) 


The measured perihelion shifts are 

^!jS ^ OBSERVED = 1+3 ± > 

OBSERVED ^ 1 > ( 15 ) 


where there may be an additional discrepancy in Mercury's perihelion shift 
due to a contribution (< 4" per century) of a possible solar quadrupole mo- 
ment. By combining equations (l4) and (15) we can eliminate the term involv- 
ing 7 and 8, and obtain a third limit on the parameters Q^, and Q^: in 
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order that the PPN perihelion shifts agree with the measured shifts within 
the experimental error, the parameters CC^, and a must satisfy 

1 13 OL x + 0.07 a 2 - 170,000 a 3 | < .2 . (16) 

We now combine the three experimentally determined limits, equations 
(8), (12), and (16) and obtain individual upper limits on the values of 
a i’ a 2’ an< ^ 3S Table 1. For comparison, Table 1 also lists 

the predicted values for CC Ot^, and Q:^ for all the metric theories which 
our groups at Caltech and Montana State have examined, and which, until 
now, were considered viable, i.e. agreed with the "classical" tests [see 
Ni ( 1972a) and Thorne, Will, and Ni (1971) for examples of "nonviable" 
theories] . 

From Table 1, we can see immediately that the stratified theories due 
to Page and Tupper, Yilmaz, Papapetrou, Ni, Coleman, and Rosen cannot be 
correct theories of gravitation — they disagree violently with experiment . 

In fact, we can be more general. We have shown (Paper i) that any Stratified 
Theory (with time-orthogonal, conformally flat space slices) which agrees 
with light-deflection and time-delay experiments must have (X^ ~ - 8, which 
is forty times larger than our experimental upper limit. Hence no "Stra- 
tified Theory" of gravity, past, present, or future, can be the correct 
theory of gravity . 

A second conclusion emerges from Table 1: the squared magnitude of the 
cosmological vector field (k) in the Vector-Metric Theory of Paper I must 
satisfy 

K 2 < 3 X 10' 2 , (17) 

in order to agree (within the experimental uncertainty) with Earth- tide data. 
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From this point of view, the results we have obtained in this paper are 
complementary to the Hughes-Drever (isotropy of inertial mass) and "ether- 
drift" experiments. Those experiments put limits on the strengths of cos- 
mological vector or tensor fields which couple to matter's nuclear or 
electromagnetic energy (Peebles 1962, Peebles and Dicke 1962, Dicke 196U), 
while our results put limits on vector or tensor fields which couple to 
matter's gravitational energy. There is a wide class of such vector-metric 
and tensor-metric theories (Paper i), and the limits on Q^, OL^, and a ' 
should be pushed as low as possible in order to put more stringent limits 
on these cosmological vector and tensor fields. A future paper in this 
series may study these theories and their predictions in detail. 

Recent research by Ni has shown that it may be possible to push the 
limits on Ot CL g, and a ^ as low as 10 by means of studies of the pulsa- 
tions of white dwarf stars. According to the PPN formalism, motion of 
pulsating white dwarfs relative to the Universe should produce instabilities, 
i.e. exponential growth or decay of the amplitudes of their pulsation. 

Results of observational studies of white dwarf pulsations, proper motion and 
stabilities may put stringent limits on such velocity- induced effects, which 
thereby put limits on CUp Q^, and a (Ni 1972b). 

A third conclusion obtained from Table 1 is that this paper's results 
do not distinguish between general relativity and scalar- tensor theories: 
these theories are not preferred- frame theories and predict no preferred- 
frame effects (for a discussion see Paper i). 

The remainder of this paper gives detailed calculations of preferred- 
frame effects. In §11 we derive the expression (eq. [3]) for the Newtonian 
Gravitational Constant as measured by a gravimeter on the Earth, and compute 
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the amplitudes of Earth-tides and variations in the Earth's rotation rate 
produced by the variations in G. We derive in §11 the perturbations on 
planetary motions produced by the additional metric terms (eq. [l]), and 
calculate the resulting secular changes in planetary orbital elements and 
the resulting periodic perturbations in the Earth-Moon distance. In §IV 
we show that other reasonable assumptions for the direction of w do not 
significantly affect our limits on Ct v a 2 , and <V Concluding remarks are 
made in §V. In an Appendix we compute the amplitude of the spherical 
expansion and contraction ("breathing") of the Earth produced by t he varia- 
tion in G. 

II. GEOPHYSICAL PREFERRED- FRAME EFFECTS 

We consider a test body (gravimeter) maintained a constant proper 

distance R from the center of the Earth by a four-acceleration F (supporting 
P 

force of the ground). Then the gravitational constant G measured by this 
gravimeter is related to F^, the radial component of this force (toward 
the center of the Earth) by Newton's law of gravitation, written in invari- 
ant language (see Will [1971] for discussion) 

F r = GM/R p 2 + F c , (18) 

where M is the mass of the Earth, and F £ is an invariant expression for the 
centrifugal acceleration due to the Earth's rotation. Will (1971) obtained 
an expression for G by evaluating equation (18) in a coordinate system at 
rest with respect to the preferred Universal rest- frame. The new version 
of the PPN formalism discussed in Paper I, allows us to evaluate equation 
(18) in the momentary rest-frame of the Earth. In this frame, the PPN 
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metric contains additional terms (eqs. [1]) which depend on w the Earth's 
velocity relative to the preferred Universal rest-frame. We wish to con- 
sider only the effects of these terms; post-Newtonian effects on G due to the 
Sun and planets and to the Earth's self-gravitational field have been dis- 
cussed elsewhere (Nordtvedt 1971; Will 1971). Other post- Newtonian effects 

-9 

produce gravimeter-measured forces smaller than 10 g (Will 1971). These 
w-dependent terms enter equation (18) only via the equation of motion of 
the gravimeter, 


, i / , _i j k i 
du /dT + r ^ u u = F , 


(19) 


where u is the gravimeter's four-velocity and t is proper time along the 
gravimeter's world line. For a gravimeter at rest on the surface of the 
Earth, equation (19) simplifies to (neglecting forces smaller than 10 ^ g). 


a a „ JL 

1 00 * 2 8 oo,a 


( 20 ) 


We note that, to our order of approximation. 


R — R — | X — X , 
p 1 ~© ~GRAVIMETER 


and use equations (l) along with the Newtonian part of g to obtain 


F r = ^ 8 


00 


(£/r) ' v[u - ft*, - a 2 - a 3 ) u - i a 2 ^ 


• ( 21 ) 


For a spherically symmetric Earth (the Earth's oblateness has significant 
Newtonian effects, but negligible post-Newtonian preferred- frame effects), 
we have 
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u = m/r , 


( 22 ) 


V = “a V r3 - 1( i / r5 > < 3 R a R p - r2 V , (23) 

where I is the spherical moment of inertia of the Earth. Equations (22) 
and (23) along with (21) give (ignoring the centrifugal acceleration): 



where 

£ r = R / R • (25) 

From equation (2*0 we obtain the gravimeter-measured value of G as given 
in equation (3). Except for the contribution due to the Earth's moment of 
inertia I, this result is in agreement with the result of Will (1971). 
Because of the Earth's orbital and rotational motion, this measured value 
of G changes (as w and § change) according to (cf. Eq. [5]) 

&/G = (| « 2 + « 3 - «i) w • v 

+ 1 a 2 [(w ‘ S r ) 2 + 2(w • e r ) (v • e r ) + (v • e r ) 2 ] , (26) 

where v is the Earth's orbital velocity, and w is the solar system's uni- 
form velocity through the Universe, which we have assumed to be 200 km/sec 
in the direction = 90°, b^ = 0^. In terms of the Geocentric Ecliptic 
coordinate system (z-axis normal to the Earth's orbit, x-axis directed 
toward the Sun at vernal equinox) the direction of w is given by A = 3b6°, 
P = 60°, and in the Geocentric Equatorial coordinate system (z-axis normal 
to the Earth's equator, x-axis same as before), it is given by a = 318°, 
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8 = ^8° (see Smart [1960] for definitions of these astronomical coordinate 
systems, and for equations to transform from one to the other). In obtain- 
ing equation (26), we have used the fact that, for the Earth, 


I - -§ m2 > (27) 

and have assumed a circular Earth orbit (v constant). 

In order to compare this variation in G with geophysical data, we must 
perform a harmonic analysis of the terms in equation (26). The frequencies 
involved will be the sidereal rotational frequency of the Earth fi, due to 
the changing e^ relative to the fixed direction of w, and its orbital 
sidereal frequency cj due to the changing direction of v relative to w, 
along with harmonics and linear combinations of these frequencies. We 
work in Geocentric Ecliptic Coordinates, and assume a circular Earth orbit, 
with the Earth at vernal equinox at t = 0. Then 


w = w cos p(cos A e + sin A e ) + sin P e I , 

~ ~x ~y ~z I ’ 


v = v(sin wt e - cos wt e ) 
~ ~x ~y 


(28) 

(29) 


For a gravimeter stationed at Earth latitude L, 


cos L cos(fit - e) 

+ I cos L sin(fit - e) cos 0 + sin L sin 0 



- cosLsin(fit - e) sin 0 


sin L cos 0 e 


(30) 


where e is related to the longitude of the gravimeter on the Earth, and 0 
is the "tilt" (23 l/2°) of the Earth relative to the Earth's orbit (ecliptic). 
Equations (28), (29), and (30) give 
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w • v = wv cos p sin(ojt -A) , 

(w • e r ) 2 = w 2 [| + §0 - sin 2 s)(| - Sin 2 l) 

+ sin 2g sin 2L cos (fit - e - O') 

12 2 1 

+ cos S cos L cos 2(fit - e - Ct) I , 

(w • e )(v • e ) = wv ~ cos p sin(uit - A) 

+ ^ - sin 2 l)| cos P sin(ut - A) + ~ sin 5 sin 0 cos cotj 
+ ^ sin 8 (l - cos 0) sin 2L sin[ (fi + to) t - e ] 


(31) 


(32) 


(Y 


- -jj- cos 5 sin © sin 2L cos[ (fi + to) t - e - CC] 

1 

- -jj- sin 5 (1 + cos 0) sin 2L sin[ (fi - to) t - e] 

- cos 5 sin 0 sin 2L cos[ (fi - to) t - e - a] 

1 2 

+ cos s(l- cos ©) cos L sin[ (2fi + to) t - 2e - a] 

1 2 ! 

- cos 5 (1 + COS 0 ) cos L sin[ (2fi - to) t - 2e - a] 


x2 2 , 

e ) = v 

~r 


3 

F 


| + |0 - sin 2 l)(| - | sin 2 ©) 

/l 2 \ 2 

l-j - sin LI sin 0 cos 2iot 


(33) 


+ sin 20 sin 2L sin(fit - e) 

- sin 0(l - cos 0) sin 2L sin[ (fi + 2to) t 

+ ^ sin 0(l + cos 0) sin 2L sin[ (fi - 2co) t 


€] 

6 ] 


(3U) 


1 ^ 



12 2 

+ |j- sin © cos L cos 2 (fit - e ) 

1 2 2 

- -g (l - cos ©) cos L cos[2(ft + to) t - 2e] 

12 2 

- — (l + cos e) cos L cos[ 2(ft - w) t - 2e] , (3^ con't.) 

where we have used both the ecliptic coordinates (A, 8 ) and the equatorial 
coordinates ( a , , 5 ) corresponding to the direction of w, in order to simplify 
the various expressions. 

Equations (31), (32), (33), and (3^) reveal four different types of 
variations in G. 

i) Semi-Diurnal Variations : These are the terms which vary with frequency 

around 2 ft: 2 ft, 2 ft + to, 2 ft - u, 2 (fi + to), 2 (ft - to); i.e. have periods around 

2 

twelve hours (to « ft) and vary with latitude according to cos L. These 
variations are completely analogous to the twelve-hour solid-Earth tides 
produced by the Sun and Moon, called "Semi-Diurnal Sectorial Waves" by 
Melchior (1966). The true gravimeter measurements for these tides are 
affected not only by the variation in G, but also by the displacement of 
the Earth' s surface relative to the center of the Earth, and by the defor- 
mation of the Earth. This variation in gravimeter readings is related to 
the variation in G by 

(Ag/g) SEMI _ DIURML = i. 18 ^/ g ) semi _ diurnal > (35) 

where the factor 1.18 is a combination of so-called "Love Numbers", which 
depend on the detailed structure of the Earth (Melchior 1966). 

ii) Diurnal Variations : These are the terms which vary with a frequency 
around ft: ft, ft + w, ft - to, ft + 2w, ft- 2to; i.e. have periods around 2k hours, 
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and vary with latitude according to sin 2L. These variations are com- 
pletely analogous to the 2k - hour "Diurnal Tesseral Waves" of the solid 
Earth (Melchior 1966), and give gravimeter readings related to the varia- 
tion in G by the same factor: 

(Ag/g) DIURNAL = 1 - 18 (^ j / g ) diurnal • (36) 

iii) Long- Period Zonal Variations : These are the variations with fre- 

2 

quencies to and 2co, and with latitude dependence (l/3 - sin L), which are 

completely analogous to the long-period tides produced by the Sun and Moon, 

called "Long-Period Zonal Waves" by Melchior (1966). These long-period 
zonal waves produces variations in the Earth's moment of inertia, which in 
turn cause variations in the rotation rate of the Earth. These rotation- 
rate variations are related to the amplitude of the zonal variations by 
(Mintz and Munk 1953; Melchior 1966) 

(Afi/fi) Z0NAL = O.Ul A zqnal , (37) 

where Atonal is related to the zonal variations in G in equations (33) and 
(3*0 by 

(Z£/g) zqnAL = a zqnal (^ ~ sm l) . (38) 

iv) Long- Period Spherical Variations : These are the variations which 
have frequency to, but no latitude dependence; they represent a yearly varia- 
tion in the strength of G, and have no counterpart in Newtonian tidal theory. 
These variations produce a purely spherical deformation of the Earth, as 
opposed to the Sectorial, Tesseral, and Zonal waves which produce purely 
quadrupole deformations. This yearly spherical "breathing" of the Earth 
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as G varies causes a variation in the Earth's moment of inertia, which in 
turn causes a variation in the rotation frequency, given by 

(^/n) spHERICAL = - (ai/i) - (1/10) (ag/g) spherical ’ 

Detailed calculations of this change in the Earth's moment of inertia 
due to the spherical variation in G are given in an Appendix. 

By combining equations (3l), (32), (33), and (34) with the expression 
for A3/G, equation (26), substituting numerical values 

1 ° 

v - 30 km/ sec , w - 200 km/ sec , © - 23 — , 

A - 346° , a - 318° , 

P - 60° , 5 - h 8 ° , (40) 

and using equations (35), (36), (37), and (39), we may compute the ampli- 
tudes of all the various components of the Earth tides (Ag/g) and of the 
variations in the Earth's rotation rate (A2/fi). These amplitudes are listed 
in Table 2. The largest predicted Earth-tide components are the sidereal 
diurnal and semi-diurnal (fi and 2n) tides. Other, smaller components 
include diurnal and semi-diurnal solar-time tides [ (ft - co) and 2 (fi- uj)], 
and tides with frequencies 2fl - w and ft - 2co. We have used the semi- 
diurnal sidereal tide ( 2 fi) to put an experimental upper limit on the value 
of a 2 (eq.[ 8 ]) rather than the larger diurnal sidereal tide (q ) , because 
agreement between Newtonian theory and observation is not as good for the 
diurnal as for the semi-diurnal tidal components, possibly because of 
diurnal effects due to heating by the Sun (Harrison et al . 1963). 

Although the zonal yearly variation in the Earth's rotation rate 
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(Table 2) is apparently larger than the spherical variation, we have ignored 

-2 

it, because the limit set on a. by the Earth tides (3 X lO - ), makes 
this effect too small to be discernible. Thus we have focused on the 
spherical variation in fi in order to set a limit (eq. [12]) on the combi- 
nation (2/3 a g + a 3 - a^) . 

Future experimental studies of Earth tides should concentrate on 
separating the various predicted components of the diurnal and semi-diurnal 
tides. Since the various components in each group are close to each other 
in frequency, a complete separation of the components would require at 
least a year of continuous gravimeter data. Such detailed studies should 
be used not only to justify our rather heuristic discussion of the limit 
set on a 2 by Earth- tides, but also, perhaps, to improve this limit. 

One further geophysical effect produced by motion through the Universe 
should be mentioned: a yearly precession or "wobble" of the Earth's axis 
of rotation which depends on the PPN parameter (X^. However, the ampli- 
tude of this wobble is too small (< 10 seconds of arc) to be discerned 
from the Newtonian Chandler wobble. 

III. ORBITAL PREFERRED-FRAME EFFECTS 

In this section, we derive a variety of effects on the orbital motions 
of planets and satellites, which may result from our motion at 200 km/sec 
through the "preferred Universal rest- frame". We consider a two-body system, 
one a "test" body, the other a massive body, with self-gravitational energy 
2 and spin-angular velocity A, moving through the Universe with velocity w. 
Using the PPN metric written in the "rest" coordinate system of the two 
bodies (thus containing the w-dependent metric terms equations [ 1 ]), we 
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derive those perturbations in the equations of motion for these bodies, 
which depend on w and on the PPN parameters a p and Q^. We then 

specialize to the solar system, and calculate (a) secular changes in the 
orbital elements of a test body (planet or satellite) orbiting the Sun, 
and (b) periodic perturbations in the relative distance between the Earth 
and the Moon. 

The perturbation 5a T , in the instantaneous acceleration of the test 
body due to the w-dependent metric terms (eq. [1]) is calculated from the 
geodesic equation: 


e, a CC,, * -rPt i j „0 i j a. 

8a T = 8(dv /dt) = ~5 r i j v v J + S T ^ ^ v v J v , 

where 5 p 1 is the contribution of gg (eq. [l]) to the Christoffel 
J K i J 

symbols. Equations (l) and (ho) yield 


m 


Mr 


^<*2 +0 5 - - I a 2 ( : ~ r -) 


Mr 1 ^ 

’ ( 2 “i 1 + “2 s) 


(H) 


where M is the mass of the central body, r is a vector from the central 
body to the test body, and v is the test body's velocity. In calculating 
§a T we have idealized the central body as a point mass: any contributions 
due to the central body's finite size (moment of inertia i) while important 
for geophysical effects at the Earth's surface (§Il), have negligible direct 
effect on the motion of orbiting test bodies. However, the massive body 
itself, although not affected by the test body (which has negligible mass) 
is perturbed by its own motion through the Universe. This perturbing 
acceleration is due to a coupling between the massive body's 


19 



self-gravitational energy, its rotation, and its velocity through the 
Universe. The test body, with negligible self-gravity, is not affected 
by this coupling. 


We compute sa^ as follows: each element of matter (with mass tn^) in 
the massive body is perturbed by the body's motion through the Universe 
according to equation (40), where now the gravitational potentials which 
appear in the Christoffel symbols are generated all other matter in the 
massive body. Using equations (l) and (4o) along with the expressions U, 
U^ and (Paper i), we get, for the i-th element of mass 


5a. = - 2 

i . / 


m . r . . 
J 


j/i r - 
r 


^ a 2 + a 3 ~ °l) ^ " I^ a i ‘ ^2 " ~ ’ ~j 


- | <V” • - 1 a 2 (s • EijAi/ - 3 a 2 ( s - r..)(v. • r^/r^ 


m. r. . 

- a, £ • wv. 

2 jA . tJ s 

i m. r, , 

+ .5. ' (a i 

JA r.. 


2a.) v. - a v. - 2a 0 w 

2 ~ j 1 ~i 2 ~ 


(1+2) 


The motion of the center of the mass of the body is given to sufficient 
accuracy by, 


\ = 2 m. a. , M = 2 m ± , 
i i 


(43) 


so the perturbation $a^ is given by 


M 8^ = 2 m t & a. 


(44) 


Combining equations (42) and (44), and simplifying the resulting expressions, 
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we get 


M 8^° = " a 3 H ap " a 2 wP ^ d n ap/ dt ^ ■ ^ a i " a 2* w ° ! ( dn / dt ) ’ (te) 


where ft^ and ft are the massive body's self-gravitational energy tensor 


and scalar respectively, given by 

a P 

, m . m . r . . r 

n = - — 2 - f . J U il 

u ap 2 . 3 ’ 

ij r. . 

ij 


1 m i m i 

a = = Q oa = " 2 * -777 ' (46) 

a iJ iJ 


and is given by 


H_ = E 

ap 

J-J 


a 6 

m. m . r . . v . 

1 -1 ij .1 
3 

r. . 
ij 


(V7) 


For a body in static equilibrium, ft^ and ft are constant, so we get, from 
equation (45), 


M 83 m“ * - °3 H CO • 


(>t8) 


We now idealize our central massive body as a nearly spherical body, 
rotating uniformly with angular velocity A. For our purposes, this is a 
reasonable model for the Sun (oblateness less than three parts in 10^) or 
for the Earth (oblateness less than one part in 10 ) . Then 


Yj = l* (Jj ' * c } ’ 


a opr 

or v . = e 




(h9) 


where x c is the position of the center of mass of the massive body. 

7 U 


Then 


e P7S A ft 

7 sa 


(50) 
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For a spherical body, we have 


‘Qp 


1 
3 


V ft J 


H 


m 


l \ « 


Thus, from equations (b8) and (52), 


S% = ' | a 3 fa/M) w J( A 


(51) 

(52) 


(53) 


The relative perturbing acceleration 5a between the massive body and the 
test body is defined by 


5a = 8a T - 5^ , 


and from equations (hi) and (53) is 

|(a 2 + a 3 - Oj) w 2 - | ctjfw • v) - | 

s ' (? “l X + a 2 ?) + 1 “ 3 ( 1 ) S 5 i 
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Mr 

~ 3 


- w 


(5h) 


We now calculate the effects of §a in the solar system. 

a) Secular Preferred-Frame Orbital Effects 
We consider a planetary orbit with the following instantaneous orbital 
elements: eccentricity e, semi-major axis a, and angle of perihelion relative 
to the equinox £j. We assume the orbit has zero inclination i relative to 
the ecliptic. 

Following the standard procedure for computing perturbations of orbital 
elements (Smart [1953], Robertson and Noonan [1968]), we resolve the per- 
turbing acceleration 5a (eq. [5h]) into a radial compontent ft, a component 
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normal to the orbital plane, and a component § normal to y and R, and 
calculate the rates of change of the orbital elements using the formulae (in 
the notation of Robertson and Noonan [1968]): 


dto 

dt 


pR 

he 


cos 


cp + 


g (p + r) 

he 


sin cp 


(55) 


If = -Mr- [ a s sln <p + f (? - r )] < (56) 

ft - tt(t + sln 1>) > (57) 


where h is the angular momentum per unit mass of the orbit, cp is the angle 
of the planet measured from perihelion, and p is the semi-latus rectum 
given by 



(58) 


We calculate the perturbations (eqs. [55], [56], and [57]) to first order 
in 5a ? retaining only secular terms, and using a Keplerian ellipse as 
unperturbed orbit, given by 


r 


p(l + e cos 



(59) 


r 2 (dcp/dt) 


= h = constant 


(60) 


We also assume the Sun's spin axis (A) is normal to the planet's orbital 
plane. For the secular changes over one orbit, our results are, to zero' th 
order in the eccentricity (e / 0) : 
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Aw 


Ae 

e 



( 61 ) 


(62) 


Aa = 0(e 2 ) 


(63) 


where w p and w^ are the components of w in the direction of the planet's 
perihelion (w p ) and in the direction at right angles to this (w^), in the 
plane of the orbit. The perturbation &a of equation (5b) can also be shown 
to produce secular changes in the inclination and angle of nodes of orbits, 
proportional to the component of w normal to the orbital plane. 

In evaluating the perihelion shift, equation (61) for Mercury and the 
Earth (cf. eqs. [lb]), we have used standard values for the orbital elements 
(Allen 1963), our adopted value for w (200 km/sec in the direction A = 3b6°, 
f3 = 60 ), and numerical values for the Sun's gravitational energy and rota- 
tional angular velocity: 

( |n!/M) 0 ax 10“ 6 , A 0 « 3 X 10~ 6 sec" 1 . 


b) Periodic Preferred-Frame Orbital Effects 
Our motion through the "preferred" Universal rest- frame may produce a 
variety of periodic perturbations in planetary orbits. Of particular interest 
are the predicted perturbations in the Earth-Moon distance, since such 
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perturbations can be studied by means of laser- ranging to the Moon. 

We idealize our Earth-Moon system as consisting of a test body (Moon) 
in orbit around a massive body (Earth) with mass M and self-gravitational 
energy The Earth is assumed to rotate uniformly with angular velocity 

A with its spin axis inclined at an angle 0 ("tilt" of the Earth) relative 
to the ecliptic. We assume the Moon's unperturbed orbit is in the plane 
of the ecliptic (the actual inclination of the orbit is only ~ 5°), and 
is given by a nearly Newtonian ellipse with a small perigee advance of 
2rte per orbit: 

r Q = p 1 1 + e cos (l - e) cpj 1 , (6*0 

h Q = r ^ qj = constant = (Mp)*'/ 2 , (65) 

where the subscript Q refers to unperturbed quantities. The perigee advance 
of per orbit is assumed to be produced by the effects of the Sun and the 
other planets and corresponds to the Moon's perigee rotation with a period 
of about 17 years, relative to the fixed stars. 

We now compute the effects of the perturbation 5a (eq. [5^]) on the 
relative distance between the two bodies. We define 

u h (l/r) = u Q + 5u = (l/r Q ) + 5u , (66) 

h = h Q + 5h . (67) 

Then to first order, the differential equations for the perturbations 5 u 
and 5 h become 

(d/dcp) 5h = (r 0 /h Q ) 2 (h Q X r Q ) * 65 , (68) 
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( 69 ) 


8u + (1 - 


e) 8u 


dcp 


■ ■ <¥ t • * - e 


dr 


° (8h) 


h o r o 


2 dcp dcp 


We then substitute equation (54) into equations (68) and (69), and 
solve these equations, using Geocentric Ecliptic Coordinates (§Il) and using 
the unperturbed expressions for r^, v^, h^ in the right-hand sides of equa- 
tions (68) and (69). We first integrate equation (68), then substitute 
our result into equation (69) and solve for su. We retain terms up to first 

order in the Moon's orbital eccentricity (e ~ .055) and to first order in 
2 

the quantity e /e ( ~ .7 for the Moon) and neglect terms which produce 
perturbations smaller than 30 cm, the current experimental limit of laser- 
ranging technology. Our result for §u is 

p S u - - |-| a^M/p) 1 / 2 (e/e) wj sin(ecp - cpj_) 

+ |a 3 ( |fl |/M)(Ap 2 /M)(e/ € ) w 2 j sin(e<p - cp 2 ) 

+ w l 2 ] cos (2ecp - 2^) 

- ^(M/p^d/e) wj sin(cp - cp x ) 

+ [-| a 3 ( |n]/M)(Ap 2 /M)(l/e) w 2 ] sin( cp - cp 2 ) 

+ [h a 2^ e / e ^ w l 2 ) cos [( 1 + e) «P - 2 9 i] 

1 1 2 1 

J2 a 2 w i cos ( 2c P “ 2 cPi) 

- a 3 (|n|/M)(Ap 2 /M) ew 2 ] sin[(2 - e) cp - cp 2 ] , (70) 
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where w^ is the projection of w onto the plane of the orbit, given by 

= w cos 3 , (71) 

and is given by 

2 2 2 1 /2 

w^ = w cos s(sin a + cos a cos £ ' 0) . (72) 

The phases cp^ and cp^ are related to the direction of w i" a manner which 
we will not quote here: it is the amplitudes of the terms in equation (70) 
which we are interested in. 

We substitute into equation (70) the standard numerical values for 
the Moon's orbital elements: 

e - 5.5 X 10' 2 , p - 3.8 X 10 5 km , e - 4.4 x 10“ 3 ; 

for the Earth's self-gravitational energy, and rotational angular velocity: 

( |n |/ m) - 5 X lO" 10 , A ~ 7 X 10 -5 sec -1 ; 

© ® 

and for the velocity w and the "tilt" of the Earth' 

1 ° 

w = 200 km/ sec , 0 = 23 — , 

3 = 60° , a = 318° , 5 = kQ° 

The resulting amplitudes of the periodic terms in equation (70) are sum- 
marized in Table 3. 

Table 3 is by no means a complete list of preferred- frame perturbations 
in the Earth-Moon distance. Higher order expansions in eccentricity may 
yield sizable perturbations with different frequencies than those listed in 
Table 3. Because the Earth orbits the Sun, the Earth's velocity w is not 
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precisely constant but various with a period of a year. This variation 
will introduce further terms in the list of perturbations, whose frequencies 
will be combinations of the Moon' s natural sidereal frequency (argument cp) 
and the Earth's orbital frequency. Further study of these terms is under- 
way. 

Our result for the term with argument 2cp (period 13.6 days, amplitude 
1+ m) is in rough agreement with the results of a heuristic calculation by 
Vinti (1971). 


IV. THE SOLAR SYSTEM'S VELOCITY THROUGH THE UNIVERSE 

Throughout this paper we have assumed that the solar system moves through 
the Universe with a velocity w equal to its nearly circular orbital velocity 
around the Galaxy (~ 200 km/sec in the direction f 11 = 90°, b 11 = 0°). A 
more realistic value for w would be the solar system's velocity through the 
cosmic microwave radiation, when it is ultimately measured with confidence. 

Current measurements of this velocity, obtained by studying the anisotropy 
in the measured temperature of the microwave radiation (caused by the Doppler 
shift) are not yet completely reliable (Conklin [ 1969]; Boughn, Fram, and 
Partridge [1971]). Those results which have been obtained, however, are in 
rough agreement with measurements of the solar system's velocity relative 
to clusters of galaxies, obtained by studies of galactic redshifts (de Vaucouleurs 
and Peters 1968), and suggest a net velocity of ~ 200 km/sec in the direction 
l 11 ~ 290°, b 11 ~ 2k° ( see Sciama [1971] for a discussion). But because of 
the experimental uncertainties, this value for our velocity relative to the 
Universe should not be given much weight at this time. 

For the sake of illustration, however, we have repeated the calculations 
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of this paper using this new value for w. The resulting limits on the 
parameters a Ct^, CL^ are 

laj < 0.1 , \a 2 \ < 2 X 10“ 2 , |« 3 | < 8 X lo -6 , (73) 

which are not significantly different than the limits given in Table 1. 

V. CONCLUSIONS 

The PPN formalism has been used to study possible solar-system effects 
due to our motion relative to a "preferred Universal rest- frame". Some of 
these predicted effects were found to be in strong disagreement with experi- 
mental data, with the result that the PPN parameters a a 2 , and CL '■ must 
satisfy the following experimental constraints: 

laj < 0.2 , |a 2 | < 3 X l(f 2 , |C^| < 2 X 10" 5 . (7k) 

These results rule out the metric theories of gravity due to Page and Tupper, 

Yilmaz, Papapetrou, Ni, Coleman, and Rosen, and put an upper limit on the 
strength of a possible cosmological vector gravitational field in a vector- 
metric theory discussed in Paper I of this series. 

The authors wish to thank the National Science Foundation and the 
Montana State University Physics Department for their support and hospitality 
during the 1971 N. S.F. Summer Workshop in Selected Topics in Theoretical 
Physics (GZ-1919-NSF), where the initial research for this paper was carried 
out. 
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APPENDIX 


VARIATION IN THE EARTH'S MOMENT OF INERTIA DUE TO 
A VARYING GRAVITATIONAL CONSTANT 

Here we derive the variation in the Earth’s moment of inertia Al/l 
caused by a variation in the Newtonian gravitational constant AG/G. We 
focus on the part of AG/ G which is independent of position on the Earth 
(eqs. [31] and [33]), and compute the amplitude of the spherical "breathing" 
of the Earth as G varies. This calculation is similar in spirit and intent 
to the work of Murphy and Dicke (196^) who studied the effects of a uniformly 
decreasing gravitational constant on the Earth. 

Since G varies so slowly (period of one year), we will assume that the 
Earth is in hydrostatic equilibrium at each moment of time, and changes only 
quasistatically. For hydrostatic equilibrium, 

_ G( t) pm(r) , \ 

dr " ' 2 * K } 

r 


where p is the pressure, p the density, G(t) the time-varying gravitational 
constant, and m(r) the mass inside radius r, given by 


m(r) 



(A2) 


We will use m(r) instead of r as independent variable, then from equation 
(Al), we have 


p(m) 



(A3) 
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where M is the mass of the Earth. Note we have used the boundary condition 
that the pressure vanish at the surface of the Earth, i.e . , 

p(M) = 0. (A h) 

As G(t) changes, the pressure distribution changes, causing a change 
I = |e in the position of each element of matter. For a given shell of 
matter, the mass inside that shell is constant by conservation of mass. 

Then if G(t) changes by AG, we get from equations (A2) and (A3), following 
each element of matter 


Am = 0, 


(A5) 



(A6) 


But the volume of each element of matter changes, and this change can be 
related to the change in pressure using the bulk modulus H (we ignore 
temperature changes; see Murphy and Dicke [196^]): 


Ap = - h (AV/V) = - h V • | = - ( K / r 2 )(r 2 |) . 

From equation (A7) we obtain an expression for % in terms of Ap: 

- r 2 

1 f r' Ap' . , 


(A7) 


i = 


_1_ f r 'A 

4 


dr' 


(A8) 


The spherical moment of inertia is given by 

M 

I* 

1=1 4 dm. 


i 


(A9) 


and the change in I caused by the displacement of each shell of matter is 
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AI 


= 2 



r | dm . 


(A10) 


Substituting equation (A8) for % into equations (A6) and (AlO), and re- 
expressing the equations in terms of r instead of m, we get 



(All) 


(A12) 


—8 

Since AG/G is small (~ 10 ), equation (A12) can be iterated to obtain 

Ap, from which AI may be computed using equation (All) . We are interested 
only in a rough estimate for AI, so we will take only the first order 
approximation for equation (A12) . Then 


AI » - 8* f 


/'•■*/'(?) 

0 J o ' 


r' 2 dr' 


(A13) 


The radial distribution of the ratio (p/h) for the Earth can be approximated 
by the function 


(P ,/h) 


(0.26 [1- (r/R) 2 ] , 

I 0 


(r/R) s 0.2 
(r/R) < 0.2 


(Alb) 


assuming an incompressible core (Allen 1963). We also use the density 
distribution 


p(r) 


p 0 



p Q = (5/2)(3M/4*R 3 ) , 


(A15) 


which gives the observed mass and spherical moment of inertia for the Earth. 


32 



Equations (A13), (AlU), and (A15) then yield 

- ( 0 . 10 ) f . 

Other Earth models which we have tried yield similar results. 


(A16) 
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TABLE 1 


EXPERIMENTAL LIMITS AND THEORETICAL PREDICTIONS 

FOR THE PPN PARAMETERS QL , a , AND a 

1 2 3 

EXPERIMENTAL LIMITS ON C^, a , , AND (X^ 

(ajl <0.2 
|a 2 j < 3 x 10" 2 
\ct 5 \ < 2 X lo -5 


THEORETICAL 

Theory and its ^ 
adjustible parameters 

PREDICTIONS 

“l 

FOR 

a x , a 2 , and 
a 2 

a 3 

General Relativity 

0 


0 

0 

Scalar-Tensor Theories (oj,a) 

0 


0 

0 

Vector-Metric Theory (k) 

0 


K 2 /(l + 1/2 K 2 ) 

0 

Stratified Theories : 

a. Page and Tupper (a, c) 

- 4(1 + 

a) 

0 

- 2(1 + a) 

b. Modified Yilmaz 

- 8 


0 

- 4 

c. Papapetrou 

- 8 


- 4 

0 

d. Lagrangian Stratified 

- 8 


0 

0 

Theory (Ni) 

e. General Stratified 

- 8 


0 

- 4 

Theory (Ni) (p, q) 

f. Coleman (p) 

- 8 


0 

- 4 

g. Rosen (A) 

- 4(1 + 

A) 

0 

0 


See Ni ( 1972a) and Will and Nordtvedt (1972) for discussion and references. 

* 

"Stratified Theories with time- orthogonal conformally flat space slices" 
(see Ni 1972a). These theories all have the property that cn^ = - 4(1 + 7 ), 
where 7 is the PPN parameter whose value must be ~ 1.0(± .1) in order to 
agree with light deflection and radar time-delay experiments. 
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TABLE 2 


GEOPHYSICAL PREFERRED-FRAME EFFECTS 


Effect 

Frequency 


Amplitude 

i) Semi-Diurnal 

2Q - 2co 

1 

-9 

x 10 a 2 

2 

cos L 

Earth-Tides (Ag/g) 

2n - u 

1 

X 10" 8 « 2 

2 . 
cos L 


2n f 

3 

Q 

x lo" a g 

2. 

cos L 


2S7 + oj 

< 

I—* 

o 

1 

CD 



2fi + 2to 

< 

lo" 9 


ii) Diurnal Earth- 

3 

OJ 

i 

a 

< 

lO" 9 

- 

Tides (Ag/g) 

n - u) 

1 

CO 

i 

o 

rH 

X 

sin 2L 


n 

7 

CM 

s 

00 

1 

o 

T— 1 

X 

sin 2L 


n + oj 

3 

_9 

x 10 a 

sin 2L 


n + 2aj 

< 

10- 9 


iii) Zonal Variations in 

U) 

8 

-9 

x 10 a 2 


Earth's Rotation 





(m/n) 

2a) 

< 

-9 

io a 2 


iv) Spherical Variation 

t 

CO 

3 

-9 

x 10 (2/3 a 2 + - a^) 

in Earth rotation 





(a a/a) 






^Used to put limits on the values of the parameters O. y a 2 , and (§l). 
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TABLE 3 


PREFERRED- FRAME PERIODIC PERTURBATIONS 
IN THE EARTH-MOON DISTANCE - A PARTIAL LIST 


Amplitude of 
distance perturbation 
in meters 

Argument of 
periodic term 

Period 

3 m 

e cp 

17 years 

8000 QL m 
o 

e cp 

17 years 

b a 2 m 

2 e cp 

8 l/2 years 

25 m 

<P 

27.3 days 

70, 000 m 

<P 

27.3 days 

40 a 2 m 

(1 + e) cp 

27. 1 days 

b a 2 rn 

2 cp 

13.6 days 

12 a : m 

(2 - e) cp 

13.7 days 
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